BALLS-IN-BOXES DUALITY FOR COALESCING RANDOM WALKS 
AND COALESCING BROWNIAN MOTIONS 
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o 
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CNJ . Abstract. We present a duality relation between two systems of coalescing random 

fH ' walks and an analogous duality relation between two systems of coalescing Brownian 

^ . motions. Our results extends previous work in the literature and we apply it to the 

study of a system of coalescing Brownian motions with Poisson immigration. 



1. Introduction 



^.. 
r^ , Consider a system of m indexed particles with locations in M that evolve as follows. 

'jrt ! Each particle moves according to an independent standard Brownian motion on M until 

two particles are at the same location. At this moment a coalescence occurs and the 
particle of higher index starts to move together with the particle of lower index. We 
say the particle with higher index is attached to the particle with lower index, which is 
>■ I still free. The particle system then continues its evolution in the same fashion. Note that 

2^ ■ indices are not essential here, the collection of locations of the particles is Markovian in its 

fT^ . own right, but it will be convenient to think of the process as taking values in R™ rather 

^ ! than subsets of M with at most m elements. For definiteness, we will further assume that 

the particles are indexed in increasing order of their initial positions: it it clear that the 

O ' dynamics preserve this ordering. Call the resulting Markov process X = (Xi, . . . , Xm)- 

r-| ■ The analogous coalescing simple random walk has many applications. One successful 

^ ! example is in voter model, which is particularly well understood because of a duality 

relation with the coalescing random walk (see, for example, |Gri79t |Lig99| ). Similarly, 
coalescing Brownian motion plays a key role in analyzing certain complex interactive 
- stochastic systems. For example, in |DEF"'"00] the coalescing Brownian motion is dual to 

rS I the Brownian stepping-stone model in the sense that it determines the joint "moments" of 

c^ I the latter. This interplay leads to further results on the Brownian stepping-stone model 

in |Zho03j . A "continuous family" of coalescing Brownian motions, usually referred to as 
the Arratia flow, serves as a fundamental example in the theory of stochastic flows. See 
|Arr79|, lHar84 l| for accounts of this topic. The Arratia flow is an example of an interesting 
noise that is not generated by Brownian motions or Poisson processes ^Tsi98, Tsi04j . More 
general "sticky" flows have recently been considered in |LJR021 ILJR031 IWar02j . 

Closed form analytic expressions for features of the joint distribution of coalescing 
Brownian motion are rarely known, but some intriguing relationships have been observed 
for stochastic systems involving coalescing Brownian motions. A self-duality relation for 
the Arratia flow is described in the Introduction of IArr79| . where the borders between 
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clusters (that is, pre-images of particles) are shown to have the same joint distribution 
as the locations of particles. A duality between a system of coalescing Brownain motions 
and a system of annihilating Brownian motions is established in |DEF"'"nnj . A dual rela- 
tionship is presented in [STWOOj between two system of Brownian motions, in which one 
system runs forward in time, the other runs backward in time, Brownian motions from the 
same system coalesce and Brownian motions from different systems reflect on each other. 
Another result along this line is obtained in in |TW97j . which involves a duality on two 
flows of Brownian motions moving at opposite directions of the time interval (—00, 00). 
Within each flow, the Brownian motions coalesce, and meanwhile each Brownian motion 
is either reflected or absorbed at depending on when it reaches 0. There is no interaction 
between the two flows. 

The distribution of X(t) is uniquely specified by knowing for each choice of 2/1 < 2/2 < 
... < Hn the joint probabilities of which "balls" Xi{t), X2{t), . . . , Xm{t) lie which of the 
"boxes" [2/1,1/2], [1/2, 2/3], ■ • ■ , [Vn-i, Vn]- That is, the distribution of X(t) is determined by 
the joint distribution of the indicators 

for 1 < i < m and 1 < j < n — 1. 

Suppose now that Y := {Yi, . . . ,¥„) is another coalescing Brownian motion. The 
distribution of Y{t) is uniquely specified by knowing for each choice of Xi < X2 < . . . < a^n 
the distribution of the indicators 

/r(t,x):=l{x,GK(t),y,+i(t)]} 

for 1 < i < m and 1 < j < n — 1. 

Thus we can think of a coalescing Brownian motion as being a set of evolving balls 
with the distribution at time t determined by how the balls fall in a fixed set of boxes, 
or we can think of a coalescing Brownian motion as giving a set of evolving boxes with 
the distribution at time t determined by how these boxes contain a fixed set of balls. We 
show that these two points of view are dual to each other in the sense that if X(0) = x 
and Y(0) = y, then for each t > the arrays of indicators {I^(t,-y)) and (/j^(t,x)) 
have the same joint distribution. We derive this duality from an analogous, essentially 
combinatorial, fact about coalescing simple random walk. 

Special cases of the above mentioned duality were proved earlier in |XZj . Instead of 
using a discrete approximation approach, the results there were directly obtained from 
coalescing Brownian motions, and, as a result, the proofs were rather lengthy. 

Moreover, we extend the Brownian motion result to a situation where the "balls" and 
the "boxes" are allowed to originate at different points in time (rather than all originating 
at time 0). This latter extended result is then used to analyse the asymptotic behaviour 
of a system of coalescing Brownian particles in which new particles arise according to a 
homogeneous space-time Poisson point process. 

The rest of this paper is arranged as follows. Section |21 contains the preparation and 
the proof of our main result on the duality between two coalescing simple random walks. 
In section |3] we generalize this dual relationship to coalescing Brownian motions starting 
from possibly different times. Some known results are re-derived. In section |3] we further 
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generalize the duality to a model with both coalescing Brownian motion and Poisson 
migration included. 



2. Coalescing random walk 

A p-simple random walk on Z is a continuous time simple random walk that makes 
jumps at unit rate, and when it makes a jump from some site it jumps to the right 
neighbour with probability p and to the left neighbour with probability 1 — p. An m- 
dimensional p-simple coalescing random walk is defined in the same way as the coalescing 
Brownian motion of the Introduction. When p = | we just call this particle system a 
simple coalescing random walk. 

Some notation is useful to keep track of the interactions among the particles in the 
coalescing system. Let Vm denote the set of interval partitions of the totality of indices 
Nm '■= {!,••• ,n}. That is, an element n of Vn is a collection vr = {Ai{7r), . . . ,Ah{n)} 
of disjoint subsets of N„j such that IJi^i(^) — ^m and a < 6 for all a G A, b G Aj, 
i < j. The sets Aii^n), . . . Ah{n) consisting of consecutive indices are the intervals of the 
partition n. The integer h is the length of tc and is denoted by /(tt). Equivalently, we can 
think of Vm as a set of equivalence relations on Nm and write i ~7r j if ^ and j belong to 
the same interval of vr G Vm- Of course, if i ~^ j, then i ~7r k ~7r j for alH < A; < j. 

Given -n G Vm, define 

ai{n) := min Ai{n) 
to be the left-hand end-point of the i^^ interval Aj(7r). Put 

Z™ := {(xi, . . . , Xm) £ ^™ '■ Xi < . . . < Xm and Xi = Xj if i ~7r j} 
and 

Z™ := {(xi, . . . , Xm) £ ^"^ '■ Xi < . . . < Xm and x, = Xj if and only if i r^^, j}. 

Note that Z™ is the disjoint union of the sets Z™, tt G Vm- 

Write X = (Xi, . . . , Xm) for the coalescing random walk. If X(t) G Z™, then the free 
particles at time t have indices ai(7r), . . . , Q;i(^)(7r) and the i^^ particle at time t is attached 
to the free particle with index 

min{j '■ 1 < j < m, j ~^ i} = max{a;fc(7r) : akiir) < i}- 

In order to write down the generator of X, we require a final piece of notation. Let 
{e^ : 1 < i < A;} be the set of coordinate vectors in Z'^; that is, e^ is the vector that has 



I 



th 



coordinate 1 and all other coordinates 0. For tt G Vm, define a map K^ : Z™ -^ Z'*^'^^ 
by 

-'^w\^) ''^TT\Xi , . . . , Xffi) . yXailTr) , ■ ■ ■ , Xq,j/ j(7r) J 

Notice that K^, is a bijection between Z™ and {x G Z''^'^^ : xi < X2 < . . . < xk^.,^)}. For 
brevity, we will sometimes write x^r for Kt^{x). 
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Write B{1/^) for the collection of all bounded functions on 17^. The generator G of X 
is the operator G : i?(Z™) — * B(Z"^) given by 

G/(x) ■.= pJ2fo K-\^^ + e?-)) + (1 - p) 5^ / o K;\^^ - eJW) 



— 1/ 



- /(tt)/ o K^^x^, / G i?(Z"), X G Z™, TT G P„. 

This expression is well-defined, because if x G Z™, then x^, x^^ + e- and x^ — e- are 
all in {x G Z'^'') : Xi < X2 < • • • < a;;(^)}. 

Note: From now on we will suppress the dependence on dimension and write e^ as e,. 

Write Z':=Z+| = {2 + |:2GZ}. An n-dimensional g-simple coalescing random 
walk on Z'" and its generator H can be defined in the obvious way. Such a process, with 
q = 1—p, will serve as the process dual to the p-simple coalescing random walk on Z"^ in 
the following way. 

Fix X G Z*" with Xi < . . . < Xm and y G Z'" with yi < ■ ■ ■ < Vn- By analogy with the 
notation introduced in the Introduction, put 

J^(t,y):=l{X,(f)G[?/„%+i]} 

and 

/^(t,x):=l{x,GK(t),y,+i(t)]} 
for 1 < i < m and 1 < j < n ~ 1. 

Theorem 2.1. Suppose in the notation above that X = {Xi, . . . , Xm) is an m- dimensional 
p-simple coalescing random walk starting at ^ = (xi, . . . ,Xm) with Xi < . . . < Xm o-nd 
Y = (Yi, . . . ,1^) is a n-dimensional (1 — p)-simple coalescing random walk starting at 
y = (1/1, ... , yn) with yi < ■ ■ ■ < ym- Then for each t > the joint distribution of the 
mx {n— 1)- dimensional random array {^ij{t, y)) coincides with that of the rax (n — 1)- 
dimensional random array (J*~(t,x)). 

Proof. For a function g : {0, l}™("^i) — > M, a vector x G Z™ with Xi < . . . < x^, and a 
vector y G Z'" with jji < ■ ■ ■ < yn, set 

9\^] y) '■— 9[^lyi,y2]\-^'i^), • • • 5 ^ly„-i,y„]{^l)i • • • > -'-[j/i,y2] l^mj) ■ ■ ■ , ^[y„-i,y„]{^m) ) ■ 

We may assume that X and Y are defined on the same probability space (r2,jF, P). We 
need to show that 

(2.1) P[^(X,;y)]=P[^(x;Y,)]. 

For X G Z™, put ^x(-) := ^(x; ■), and for y' G Z'", put ^y(-) := g{-',y)- In order to 
establish (J2.H) . it suffices by a standard argument (cf. Section 4.4 in |EK86j ) to show that 

(2.2) G{-gy){^) = H{-g^){y) 

(recall that G and H are the generators of X and Y, respectively). 

Fix X G Z^ and y G Z'^ for some n G Pm and zu G Vn- Set 

1 
/+ := {i : 1 < i < /(tt), Xa,(vr) + - = ya,{m) for some 1 < j < l{tv)} 
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and 

r ■= {i:l<i< /(tt), Xa,(n) - - = ?/a,M for some 1 < j < l{zu)}. 

Similarly, put 

1 
J := {i : 1 < i < ^(w), 2/qjM - 2 = ^".W for some I <t < /(tt)} 

and 

-^^ := {j : 1 < j < ^(ro), ?/a,M + 2 ^ ^"»W ^°^ ^°"^^ 1 < ^ < ^(tt)}. 
Recall that Xai{TT) < • • • < aj^, (^) and yai(-uu) < • • • < Uau^^ivj)- Therefore, for each 
i E I~^ there is a unique j G J~ such that Xai{TT) + | = Vajizu) and wzce versa. Fix such a 
pair («,i). Observe that 

x':=x+ J2 e^ = K-^(x^ + e,) 

fceAi{7r) 

fceAj(ro) 
Writing 1(5) (■) for the indicator function of a set B, we are going to verify that 

(2.3) {li[y,,,yr+im4)) = ili[y'^„y'^,^,])ixe)) 

by considering all the possible scenarios. 
Given any i' G ^^(vr) we have: 

• for j' = aj{w) — 1, 

!([%•'> %'+i]) (4) = ^i\yj',yj'+i]){xi> + 1) 

= 

= i(bi',%-'+i - ^])M = ^i[y'j',y'j'+i])M, 

• for <yj{^) < j' < raaxAjlzu), 

l([%-M/i'+i])(4) = !([%•', %•'+!]) (a^i' + 1) 
= 

= !([%•' - l,%'+i - l])(a;i') = l([?/^v,?/^v^_J)(xi/), 

• for j' = max Aj(tu), 

!([%•', %'+i]) (4) = l(bi',%'+i])(a;i' + 1) 
= 1 

= !([%■' - ^,yj'+i])M = l{[y'y,yj>^^]){xi,), 

• and for / < aj{zu) — 1 or j' > m.axAj{zu), 

!([%•'' %•'+!]) (4) = l(biM/i'+i])(4) = l(bi', %•'+!] )(a;i' + 1) 
= 

= myj',yj'+i])M = H[y'j',y'j'+i])M- 

Moreover, given any i' ^ Ai{n), we have Xii ^ Xa^ij,)- Hence 
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• for i' = aji^Uj) — 1, 

!([%■'> %-'+i]) (4) = myj',yj'+i])M = l ([%•',%■'+! - l])ixi>) = l{[yj,,yj,^^]){xi'), 

• for j' = raaxAj{w), 

!([%■', %•'+!]) (4) = H[yj',yj'+i])M = !([%•' - '^,yj'+i])M = my'j',yj'+i])M, 

• for aj{uj) < j' < max Aj{zu), 

l{[yj>,yj>+i]){x'i,) = l{[yj,,yj>+i]){xi') = l([%v - l,yj,+i - l]){xi>) = l{[y'y,y'y^^]){xi,), 

• and for / < djizu) — 1 or j' > maxAj^zu), 

i(bi',2/i'+i])(4) = '^i[yj',yj'+i\)M = ^i[y'j',y'j'+i\)M- 

Combining the above observations yield ()2.3|) . 
Therefore, 

9y o K^i^n + ei) =g^o K^\y^ - ej). 
Furthermore, it is easy to see for i' ^ /^ that 

9y ° K^i^TT + Gi') =gyO K^^M 
and for / ^ J^ that 

Similarly, for any i E I~ there exists a unique j G J^ such that Xai(TT) — | = yajim) and 
vice versa. For such a pair {i,j) we have 

9y ° ^^^(xtt - ei) = ^x o K-\y^ + e^). 

Furthermore, we see for i' ^ I~ that 

9y o K~^{^7T - ei') =gyO K^^M 

and for / ^ J+ that 

^x o K^^iYv. + Gf) =g^o K^^{y^). 
Lastly, note that 

9y o K\^n) = ^(x; y) = ^x o K^\yz.) 
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and so 

Ktt) 



+ (1 - P) E (^y ° ^^'("^ - ^^) - 9y o K\^^)) 

i=l 

l{zu) 

-{l-p)J2 (^x o K^\y^ + e,) - -g^ o ^^^(y^)) 
i=i 

+ {l-p)J29y° K-\^^ - e.) - (1 - p) ^ ^x o K-\y^ + e,) 
is/- jeJ+ 

= 0, 

as required. D 

Remark 2.2. One can see from the proof that it is crucial that the random walks make 
only nearest neighbor jumps. 

3. Coalescing Brownian motion 

In this section we will show that the duality in Theorem 12.11 also holds when the coa- 
lescing random walks are replaced by coalescing Brownian motions. Coalescing Brownian 
motion can be defined similarly to coalescing random walk. This duality between two co- 
alescing Brownian motions follows if one can show the unsurprising fact that a coalescing 
random walk scaled in time and space in the usual way converges weakly to a coalescing 
Brownian motion. 

Proposition 3.1. The conclusion of Theorem \2. 1\ holds when the coalescing random walks 
are replaced by coalescing Brownian motions in the definition of {I^(t,y)) and (/j^(t,x)). 

We omit the the proof, but remark that a particularly straightforward martingale argu- 
ment proof of the convergence of coalescing random walk to coalescing Brownian motion 
can be given using the following result that parallels Levy's celebrated martingale char- 
acterization of Brownian motion (and is a fairly simple consequence of that result). We 
also omit the proof of this theorem. 

Theorem 3.2. Let X be an m-dimensional continuous process with X(0) = x, where 
Xi < . . . < Xm, O'lT'd let JF-^ denote the filtration generated by X. Then the following are 
equivalent. 
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(i) The process X. is a coalescing Brownian motion. 

(a) For each 1 < i < m, the process Xi is a Brownian motion with respect to T^ , and 
for each pair 1 < i < j < n, the process -j^{Xj — Xi) is a Brownian motion stopped 

at with respect to J-'^. 
(Hi) The process X is a continuous martingale with quadratic variation {Xi,Xj)t = t — 
Tij A t, where Tij := inf{s > : Xi{s) = Xj{s)}, I < i < j < n. 

For a function g : {0, Ij^-C"^!) ^ M, a vector x G Z"* with Xi < . . . < Xm, and a vector 
y G Z'" with yi < ■ ■ ■ < i/n, set 

S'l^i yj '■— dK^lyi,^]^-^^)! • • • 5 ^lyn^i,yn]\^l)i ■ • • ; -'-[2/1,2/2] l-^mj) • • • j ^[yn~i,yn]\^m) ) ■ 

Proposition 13.11 says that for any x and y 

P[^(X,;y)]=P[^(x;Y,)]. 

By choosing the right function g we can recover some known duahties. For example, given 
n = {Ai, ..., Ah) e Vn and yi < ... < y2h, put 

h 

Mx;y) = n n My2,-uy.A^^)^ x G W\ y G M'". 
i=i ieAj 

Then Proposition 13. II implies that 

p I n fi {x,it) G b2,-i,z/2,]} > = p i n n i^^ ^ [y2,^iit),Y,,it)]} 

which gives the duahty in Theorem 1.1 of |XZj . If we choose 

n / m \ 

^(x; y) = n ^ " n^^ ~ My2j-i,y2d''i)) , X G R", 2/1 < . . . < 2/2,r 
i=i \ i=i / 

then 

j=i i=i J U=i j=i 

Therefore, Proposition 3.7 in j XZj follows readily. 

The duality Proposition 13.11 can be generalized to one involving coalescing Brownian 
motion starting from different times. In order to state this result, it will be convenient to 
think of coalescing Brownian motion a httle differently from what we have done so far. As 
we have defined it, the coalescing Brownian motion X takes values in the space {x G M™" : 
Xi < X2 < ■ . . < Xm}- It will be more convenient to work with a related process for which 
we don't impose this condition. Given an arbitrary x G M™, let a be any permutation 
of the indices {1,2,..., m} such that Xo-(i) < Xa-(2) < • • • X(j(^rn)- Let X be an M"^- valued 
process that has the same distribution as the process (Xo-i(i), Xo-i(2), . . . , Xo-i(m)), where 

X(0) = (xa-(i) , 2;o-(2) , • • • ,X(j(rn))- It is uot difficult to see that X is a time-homogeneous 
strong Markov process. The following result is obvious. 
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Corollary 3.3. The duality in Proposition VJ. 1\ holds when the ordered coalescing Brown- 
ian motion X is replaced by the unordered coalescing Brownian motion X. 

Given {{si,xi), . . . ,{sm,Xm)) G (M+ x I^)™ with < si < S2 < . . . < Sm, define a 
process X taking values in {e} U |J^-^(M+ x M)'^, where e is the null vector, as follows. 
Let < o"! < . . . < 0"^ denote the distinct elements of (si,...,Sm) written in order. 
For t G [0,cri[, X(t) = e. For t G [cri,cr2[, X evolves as X(- — ai) under the initial 
condition X(0) = (xj : Sj = ai). Inductively, if X(t) has been defined on [0,a/i[, then 
for t G [cT/i, cr/i_|_i[ (with the convention cr^+i = oo), X evolves conditionally independently 
of {X(ii) : u G [0,cr/i[} given X(w— ) as X(- — ah) under the initial condition X(0) = 
X(ti— ) U {xi : Si = ah) (where U denotes the operation of appending one vector to the 
end of the other) . 

The following result is a straightforward consequence of Proposition 13.11 and repeated 
applications of the Markov property. 

Proposition 3.4. Let ((si, Xi), . . . , {sm,Xm)) and X be as above, and letY = {Yi, . . . , Yn) 
be a coalescing Brownian motion starting at y = {yi, . . . , yn), with yi < . . . < y-m- Then, 
for t > maxj Si, the m x [n — 1)- dimensional random array 



has the same distribution as 



l[%,%-+i](^iW) 



'^[Y,{t-s,),Yj + i{t-s,)]iXi 



Given two functions f,g -.R^ ~^R with /(t) < g{t) for all t, let D^{f,g) C [0,t] x R 
denote the region sandwiched between the graphs of / and g up to time t. That is, 

Dr{f,g) := {{s,y) ■.0<s< t,f{s) <y< g{s)}. 

Let D;r{f,g) := {{t - s,y) : {s,y) G Dt{f,g)} be the region D^{f,g) time-reversed at 
time t. The conclusion of in Proposition 13.41 is that the random array 



My,,yj+i]iXi{t)) ^ 

has the same distribution as the random array 

(l{(5„x,)GDr(>;-,^,+i)})- 

4. Coalescing Brownian motion with Poisson migration 

In this section we are going to study a particle system which can be described intuitively 
as follows. Given a time-space Poisson random measure n"*" on IR+ x M with intensity 
measure A x Lebesgue. Particles appear at the atoms of n"*". Once a particle appears, it 
starts to move. The existing particles execute coalescing Brownian motion with possibly 
different initial times. Define a set- valued process S by taking St to be the set of locations 
of those particles at time t > 0. 

The easiest way to define S formally is via the coalescing Brownian fiow (j) of Arratia 
|Arr79j . Here (j){s,t,x) for s,t,x G M with s < t is a collection of random variables with 
the properties 
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• the random map (s, t, x) i— > (j){s, t, x) is jointly measurable, 

• for each s and x, the map t \-^ (p{s,t,x), t > s, is continuous, 

• for each s and t with s < t, the map x ^-^ (p{s,t,x) is non-decreasing and right- 
continuous, 

• for s < t < M, 0(t, u, ■) o 0(s, t, •) = 0(s, u, ■), 

• for M G M, (s, t, x) 1-^ 0(s + -u, t + M, x) has the same distribution as cp, 

• for xi, . . . , Xm G M the process (0(0, t, Xi), . . . , 0(0, t, Xm))t>o has the same distri- 
bution as X started at (xi, . . . , Xm)- 

We then set 

St = {(f){s,t,x) : (s,x) G n+, < s < t}, 

where we use the short-hand notation (s,x) G 11+ to mean that (s,x) is an atom of n+. 
For any 6 > we have that almost surely 

S^n [-b,b] = {0(s,t,x) G [-6,6] : (s,x) G n+, < s < t,x G [-a, a]} 

for all a > sufficiently large: in particular, St is almost surely a discrete set and we 
identify St interchangeably with the simple point process obtained by placing a unit mass 
on each point. Using this observation, conditioning on 11+ , by Proposition 13 .41 and taking 
limits, we get the following result which characterizes the avoidance function and hence 
the distribution of St (see Theorem 3.3 of |Kal76j ) . 

Proposition 4.1. Given yi < . . . < y2n, let Y be a coalescing Brownian motion starting 
from {yi,. ..,|/2n)- Then 

r I St n\J ]y2j-i, y2j] = i = P exp Ua f^ y Y2,is) - Y2,-i{s)ds 

We can re-phrase Pr op osition 14 . 1 1 as follows. For fixed ao G M, the function 

^^ lA""(0(Oi->^)'^(O>-,ao))| = / 0(0,5,6) -0(0, s,ao)rfs, 6 > Oq, 

Jo 

where | • | denotes Lebesgue measure in the plane, is non-negative, non-decreasing, and 
right-continuous. It follows that there is a unique random Radon measure Mt on M such 
that 

Mt(]a,6])= / 0(O,s,6)-0(O,s,a)ds, 6 > a. 

Proposition 14.11 then says that St is the simple point process obtained by placing a unit 
mass at each atom of the Cox process with the random intensity measure \Mt\ that is, 
conditional on Mt = m, St is distributed as the random measure which places a unit mass 
at each atom of a Poisson process with intensity measure Xm. Note that Mt has atoms, 
and so the resulting Cox process will not be a simple point process; that is, it can have 
atoms with mass greater than one. Consequently, St is not a Cox process. 
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There is a unique random Radon measure M^^ on M such that 

poo 

M^{]a,b])= (/)(O,s,6)-0(O,s,a)rfs 

Jo 
= \imMt(\a,b]) 

t— >oo 

= sup M((]a, 6]), b>a, 
t>o 

(the finiteness of Moo(]a, b]) is assured by the continuity of s i-^ 0(0, s, a) and s \-^ 0(0, s, b) 
and the fact that 0(0, s, a) = 0(0, s, 6) for all s sufficiently large). Hence St converges in 
distribution as t ^ oo to the simple point process obtained by placing a unit mass at 
each atom of the Cox process with the random intensity measure XM^o- As with St, the 
point process 5*00 is not a Cox process. 

We can give an almost sure construction of S^o as follows. Consider a Poisson random 
measure n~ on M_ x M with intensity measure A x Lebesgue. Then St has the same 
distribution as 

{(p{s,0,x) : {s,x) eW, -t<s< 0}, 
and so 5*00 has the same distribution as 

{0(s, o,x) : (s,x) e n^, -oo < s < 0}. 

We can do some explicit computations for 6*00. In what follows, let Ai denote the Airy 
function - see |AS72j for its definition and related properties. 

Proposition 4.2. For a < b, 

Ai(x^b-a)) 

and 

P[#5o.n]a,6]] = (3A)i^(6-a). 

Proof. Let (1^1,12) be a two-dimensional coalescing Brownian motion starting at (a, 6). 
Then -i^{Y2 — Yi) is a Brownian motion stopped at 0. By Theorem 1 in |Lef89j . Theorem 
1 in |Lac93j . or Proposition 5.14 in |GS79j . we have 

F{S^n]a, 6] = 0} = P [exp(-AMoo(]a, b]))] 

P exp( -A [ ¥2(3) - Yi{s) ds 



Ai(x^b~a] 
Ai(Or" 
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Note that 



Thus, 



hm 

rf-cj.0 



P{5oon]c,(i]^0} d 



Ai I X^x 



d — c 



x=0 



-A3 



P[#^oon]a,b]]= hm Vp<^5oon 

n— +00 ' ^ 

i=l "^ 



(ia; Ai(0) 

Ai^(O) 
Ai(0) 



(i — l)(b — a) i(b — a) 
a+ -,a + - 



n 



n 



+ 



D 
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